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The three-state antiferromagnetic Potts model on the simple cubic lattice is investigated using 
the cluster variation method in the cube and the star-cube approximations. The broken-sublattice- 
symmetry phase is found to be stable in the whole low-temperature region, contrary to previous 
results obtained using a modified cluster variation method. The tiny free energy difference between 
the broken-sublattice-symmetry and the permutationally-symmetric-sublattices phases is calculated 
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' in the two approximations and turns out to be smaller in the (more accurate) star-cube approxima- 

, tion than in the cube one. 

PACS numbers: 05.50.+q, 64.60. Cn 

The three-state antiferromagnetic Potts model on the simple cubic lattice has been intensively studied in the last 
O ' two decades, because of its unusual and rather obscure low-temperature and critical properties [p]-^6|. The model 
, hamiltonian is very simple and can be written as 



X 



n = jJ2s(s i ,s j ), (i) 



^— > 

where J > is the interaction strength, the summation is over all nearest neighbour (NN) pairs, 8 denotes Kronecker's 
d • delta, and the variables Si can take on three different values, say {1, 2, 3}. 

It is now commonly accepted that the model exhibits a continuous phase transition at temperature T c = 1.23 (in 
units of J/ksi with Boltzmann's constant), belonging to the universality class of the three-dimensional XY model, 
as suggested by the value of the critical exponents recently calculated by means of extensive Monte Carlo simulations 
Ji~3"| and coherent-anomaly method Q|. The nature of the ordered phase, however, is not yet as clear as the critical 
O ' behaviour. 

In the past some authors suggested that the ordered phase should be one with algebraically decaying correlations 
' like a Kosterlitz-Thouless phase fipl[7|||], or an incompletely ordered phase [^0|. Other studies, however, indicated 
. the existence of a long-range ordered phase called broken-sublattice-symmetry (BSS) phase (2]-|4|||. In this phase 
the simple cubic lattice is split into two interpenetrating sublattices, say A and B, such that a site in A has all its 
. nearest neighbours belonging to B and viceversci. In terms of the site expectations Pk,a ~- 

(5(si,k)),i € a, with 

k = 1, 2, 3 and a — A, B, the BSS phase can be roughly described, at very low temperatures, by pi.a — P Pi.b = 
and pk,A — 0,Pk,B — 1/2 for k = 2,3 (it can be easily verified that the phase is sixfold degenerate). To be more 
\^ ' precise, even at very low temperature, there is a small probability to find Si = 2 or 3 in sublattice A, too, because 
0^ ■ this causes no increase in energy when the neighboring sites are all in the state 3 or 2, respectively. This implies that 
the long range order does not saturate even at zero temperature, and that the zero temperature entropy per site is 
slightly larger than the value \ In 2 ~ 0.346574 one would predict naively. 

Furthermore, in a recent investigation of the Blume- Emery-Griffiths model by Rosengren and Lapinskas JHJ, based 
. on a modified cluster variation method, a new phase was found in an intermediate temperature region (0.78 < T < T c ) 
between the BSS and the disordered phase. This phase, of degeneracy twelve, is characterized by the relation 
Pi,A — Pi,b > P3.A = Pz.b > Pi, a = Pi,b between the site expectations, and has been called permutationally- 
' symmetric-sublattices (PSS) phase by the authors. The result is qualitatively similar to that obtained in the Bethe 
. approximation (which, of course, predicts higher transition temperatures), and have received a partial confirmation 
by a Monte Carlo investigation by the same authors and Kundrotas jl6| . On the other hand, a different Monte Carlo 
investigation by Kolesik and Suzuki |l5[ suggests that an intermediate phase might exist but cannot be the PSS phase 
predicted by Rosengren and Lapinskas. Kolesik and Suzuki have computed the free energies of the BSS and PSS 
phases, finding indications that the BSS free energy is lower in the low temperature region, up to a temperature which 
is higher (close below 1.0) than the transition temperature proposed by Rosengren and Lapinskas. For even higher 
temperatures, the free energy difference between the BSS and PSS phases is indistinguishable from zero, although the 
authors conclude that they cannot definitely rule out the possibility that it vanishes only at the critical point. 

In the present work we examine carefully the issue of the type of ordered phase(s) by means of two high order 



approximations of the cluster variation method (CVM) 17-19]. The CVM is a simple and powerful variational 



method for Ising-like models, based on a cumulant expansion of the variational principle of statistical mechanics. The 
full (and generally unaffordable) variational principle is based on the minimization of the functional 



1 



F[pa] = Tr(p A W + k B Tp A lnp A ), 



(2) 



where A is the lattice, p A the corresponding trial density matrix and Tr stands for trace. In the CVM one chooses a 
set of maximal clusters to take into account and considers the approximate variational functional 



F[{p a ,aeM}]= Tr (P a T~i a 4- k B Ta a p a In p a ) , 



(3) 



a£M 



where M is the set of the maximal clusters and all their subclusters, p a the reduced trial density matrix for the cluster 
a, Ti. a the hamiltonian contribution associated to the cluster a (in the present case TL a = if a is not a n.n. pair) and 
the a Q 's are constant, geometry-dependent coefficients that can be easily obtained by solving a suitable set of linear 
equations The approximate functional F has to be minimized with respect to the cluster density matrices 

(which, for a classical model, are diagonal), with the constraints (normalization and compatibility, respectively) 



Trp Q = 1, a £ M 



and 



Tr a \[3Pa = P/3, C a G M. 



(4) 



In this work we shall use two different approximations in the CVM scheme, the cube and the star-cube ones. 
The cube approximation is obtained by selecting as maximal clusters the cubic cells of the lattice, and the 
corresponding functional (free energy per site) is 



f[p s ] = 3Tr{ P2 n 2 ) + k B T 



Tr(/9 8 lnp 8 ) - 3Tr(/9 4 lnp 4 ) + 3Tr(p 2 lnj0 2 ) - ^Tr(pi^ lnpi^) - ^r(p 1B lnp 1B ) 



(5) 



where pg, pi, p 2 and p\ a denote the density matrices for the cube, square, n.n. pair and site (with a sublattice index 
a = A,B) clusters, respectively. / can be regarded as a functional of ps only, since the compatibility constraints can 
be easily solved by defining the other density matrices as partial traces of p%. In principle, ps has 3 8 = 6561 different 
diagonal elements, but symmetry considerations reduce this number to 495, which (apart from the normalization 
constraint, which can be easily dealt with) is the final number of independent variables for this approximation. 
Despite this large number of variables, the local minima corresponding to the various phases are easily obtained by 
means of the so-called natural iteration method (NIM) pp| . 

The star-cube approximation |2l[] goes a step further, including in the set of maximal clusters also the "stars" 
formed by a site and its six nearest neighbours. This choice seems particularly useful here, since in view of the above 
remarks on the very low temperature properties of the BSS phase it should be particularly important to consider 
explicitly the local environment of a site. The resulting approximate functional has the form 



f[ps, p 7A , Pte] = 3Tr((0 2 H 2 ) + k B T 



Tr(/9 8 hips) + ^Ti(p 7A \np 7A ) + ^Ti(p 7B \np 7B ) - 3Tr(p 4 lnp 4 ) 



-4Tt(p4> a ^P4'a) - TTr^s lnp4's) + 6Tr(p 3A lnp 3A ) + 6Tr(p 3B lnp 3B ) - 3Tr(,02 hip 2 ) 



(0) 



where p 7a , pi' a and p 3a are the reduced density matrices (all with a sublattice index a referring to the central site) 
for the star cluster and its four- and three-site subclusters obtained by taking the central site and three or two of its 
nearest neighbors forming a solid angle or a plane angle, respectively. The functional now depends on /?§, Pta and 
p 7B (subcluster matrices being defined as partial traces, as above), but the diagonal elements of these matrices are 
not all independent because of the compatibility constraints. Taking into account the lattice symmetries one finds 
495 variables for the cube matrix, 168 variables for each star matrix, and two sets of 30 linear constraints. A more 
detailed account of our calculations will be given elsewhere Q . 

It is important to notice here that the calculation which suggested the stability of the PSS phase 0,|l^l was based 
on a modification of the cluster variation method which, although considering many-point clusters (up to the cube 
cluster, like our cube approximation), does not take into account properly the n-point correlations with n > 2. In the 
modified CVM, in fact, the density matrices are written as p a = exp(— Ha S ^), where Ha S ^ is an effective hamiltonian 
which contains only one-site and two-site terms. This results in a reduced number of variational parameters and in 
a violation of the compatibility constraints in Eq. |] (for a full description of the method see Ref. |^3| ; the authors 
of Refs. 0,^2) do not report details of the method used, leaving them to Ref. |24|], which in turn refers to p3|). 
The results for the ordered phases are then qualitatively similar to those obtained using the much simpler two-point 
(Bethe) approximation. The present work can then be regarded as an improvement of that described in jllj for two 
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reasons: first, the CVM is used in its full form, i.e. without the restriction inherent to the modified version used in 
p[ ; second, with the star-cube approximation we go one step further in the cumulant expansion of the variational 
principle. 

Let us now turn to a brief description of our results. As far as very low temperatures are concerned, results from 
the cube approximation are in perfect agreement with the modified CVM at the cube level In particular we 

have piA = 0.939067, p 2 A = P3A = 0.030467 and p 1B = 0.000117, p 2B = Psb = 0.499942, while the entropy per 
site is s — 0.366928. The star-cube approximation introduces only very small corrections, giving p\A — 0.939041, 
P2A = P3A = 0.030480 and p 1B = 0.000055, p 2B = Pzb = 0.499972, with an entropy per site is s = 0.366941. 

Also the critical temperature does not change much with respect to the modified CVM. In the cube approximation 
we obtain T c = 1.268, which is just 0.5 % higher than the modified CVM results, a difference of the same sign and order 
of magnitude of the one found in the simple Ising case p3fl . In the star-cube approximation the critical temperature 
is slightly lower, T c ~ 1.263. These results are to be compared with the most recent Monte Carlo estimate T c — 1.23 
& 

These estimates can be further improved by employing the recently introduced cluster- variation-Pade-approximants 
method (CVPAM) (2^j2^|, from which one can try to extract also estimates for the critical exponents. Using the star- 
cube results for the susceptibility for values of w = tanh(l/T) up to w max = 0.48 we have obtained (with simple 
Dlog Pade approximants) a critical temperature T c ~ 1.24 and a critical exponent 7 ranging from 1.30 to 1.33, in 
agreement with the best recent estimates yielding 7 ~ 1.31 |l3). These CVPAM results are still affected by relatively 
large sistematic errors, and an accurate CVPAM analysis would require both using larger maximal clusters and 
including corrections to scaling [ p6| . Improving in this direction is beyond the scope of the present paper, but our 
calculation of 7 can at least be regarded as a check of the accuracy of the CVM results. 

The central result of our paper, qualitatively different from the modified CVM, is however that the free energy of 
the BSS phase is always lower (although slightly) than that of the PSS phase. In Fig. |l| we have plotted the free 
energy difference Sf vs. temperature as given by the two approximations. It is important to notice that the numerical 
errors involved in the estimation of the free energy are several order of magnitudes smaller than Sf. From Fig. |l| one 
can see that Sf is very small in the cube approximation, and becomes even smaller in the more accurate star-cube 
approximation. In both approximations Sf vanishes only at the critical point, but in the star-cube approximation it 
becomes negligibly small at a temperature close below 1.0, like in the Monte Carlo simulations by Kolesik and Suzuki 
p5[ . Furthermore, the numerical values of Sf are of the same order of magnitude of those reported by these authors. 

Although it is difficult to judge on the basis of our results which is the exact ordered phase, some conclusions can 
certainly be drawn. First of all, the results found by Rosengren and Lapinskas with the modified CVM (TTJ are ruled 
out and attributed to the neglect of high-order correlations (remember that a simple Bethe approximation yields the 
same qualitative results) . This means also that the phase diagram they proposed for the Blume-Emery-Griffiths model 
p2| is incorrect and must be reconsidered. On the other hand, the present CVM predictions are compatible with the 
Monte Carlo results by Kolesik and Suzuki fH^ . The Monte Carlo simulations tends to favor the possibility that the 
free energy difference Sf vanishes below the critical point, giving rise to a new intermediate phase. The CVM results 
seems to support this possibility, the free energy difference Sf being strongly depressed, while going from the cube to 
the star-cube approximation, in the same temperature region. In order to make this explicit, we have plotted in Fig. 
H the ratio between the two values of Sf reported in Fig. |]. It can be seen that this ratio increases very quickly for 
temperatures between 0.9 and 1.0, and this could be a signal that going to larger and larger maximal clusters the free 
energy difference tends to zero. Both the Monte Carlo simulations and the cluster variation method, however, leave 
open also the possibility that the free energy difference vanishes only at the critical point, and further investigations 
are therefore welcome. 
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FIG. 1. Free energy difference between the PSS and the BSS phases, as given by the cube (dashed line) and the star-cube 
approximation (solid line) 



FIG. 2. Ratio of the free energy differences given by the cube and the star-cube approximations 
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